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00 ■ Applications of perturbative QCD to deeply virtual Compton scatter- 

ed ' ing and hard exclusive electroproduction processes require a generalization 

of usual parton distributions for the case when long-distance information is 
accumulated in nonforward matrix elements of quark and gluon light-cone 
operators. We describe two types of nonperturbative functions parametriz- 
ing such matrix elements: double distributions F(x,y;t) and skewed dis- 
i-S^ ' tribution functions J-^(X; t), discuss their properties, and basic uses in the 

QCD description of hard exclusive processes. 
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■ 1. Introduction 

The standard feature of applications of perturbative QCD to hard pro- 
cesses is the introduction of phenomenological functions accumulating in- 
formation about nonperturbative long-distance dynamics. The well-known 
examples are the parton distribution functions f p /H( x ) [1] used in pertur- 
bative QCD approaches to hard inclusive processes and distribution ampli- 
tudes (f n (x), ipiy(xi, X2, X3), which naturally emerge in the asymptotic QCD 
analyses of hard exclusive processes [2, 3, 4, 5, 6, 7]. 

The cases of deeply virtual Compton scattering (DVCS) and hard exclu- 
sive electroproduction processes [8, 9, 10, 11, 12, 13] involve nonforward ma- 
trix elements (p'\ . . . \p). The important feature (noticed long ago [14, 15]) 
is that kinematics of hard elastic electroproduction processes (DVCS can 
be also treated as one of them) requires the presence of the longitudinal 
component in the momentum transfer r = p — p' from the initial hadron 
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to the final: r + = Cp + - For DVCS and /9-electroproduction in the region 
Q 2 S> \t\,mjj, the longitudinal momentum asymmetry (or "skewedness" ) 
parameter £ coincides with the Bjorken variable XBj = Q 2 /2(pq) associ- 
ated with the virtual photon momentum q [16]. This means that the non- 
perturbative matrix element (p'\ . . . \p) is nonsymmetric (skewed), and the 
distributions which appear in the hard elastic electroproduction amplitudes 
differ from those studied in inclusive processes. In the latter case, one has a 
symmetric situation when the same momentum p appears in both brackets 
of the hadronic matrix element (p\ . . .\p). 

To parametrize nonforward matrix elements (p — r\ O(0, z) \ p) \ z 2 =0 of 
quark and gluon light-cone operators one can use two basic types of nonper- 
turbative functions. The double distributions (DDs) F(x,y;t) [9, 11, 17, 18] 
specify the Sudakov light-cone "plus" fractions xp + and yr + of the initial 
hadron momentum p and the momentum transfer r carried by the initial 
parton. The other possibility is to treat the proportionality coefficient £ as 
an independent parameter and introduce an alternative description in terms 
of the nonforward parton distributions (NFPDs) ^(X; t) with X = x + y( 
being the total fraction of the initial hadron momentum taken by the ini- 
tial parton. The shape of NFPDs explicitly depends on the parameter £ 
characterizing the skewedness of the relevant nonforward matrix element. 
This parametrization of nonforward matrix elements by J-^(X;t) is similar 
to that proposed by X. Ji [8] who introduced off-forward parton distribu- 
tions (OFPDs) H(x,^; t) in which the parton momenta and the skewedness 
parameter £ = r + /2P + are measured in units of the average hadron momen- 
tum P = (p+p')/2. OFPDs and NFPDs [11, 12] can be treated as particular 
forms of skewed parton distributions (SPDs). One can also introduce the 
version of DDs ("q-DDs" [18]) in which the active parton momentum is 
written in terms of symmetric variables k = xP + (1 + a)r/2. 

The basics of the PQCD approaches incorporating skewed parton distri- 
butions were formulated in refs.[8, 9, 10, 11]. A detailed analysis of PQCD 
factorization for hard meson electroproduction processes was given in Ref. 
[12]. Our goal in the present lectures is to give a description of the approach 
outlined in our papers [9, 10, 11, 17, 18]. 

2. Double distributions and their symmetries 

In the pQCD factorization treatment of hard electroproduction pro- 
cesses, the nonperturbative information is accumulated in the nonforward 
matrix elements (p—r \ O(0, z) \ p) of light cone operators O(0, z). For z 2 = 
the matrix elements depend on the relative coordinate z through two Lorentz 
invariant variables (pz) and (rz). In the forward case, when r = 0, one ob- 
tains the usual quark helicity-averaged densities by Fourier transforming 
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the relevant matrix element with respect to (pz) 

(p, s' | $ a (0)zE{0, z; A)4> a {z) \p,s)\ z 2 =0 (1) 

-^.) J f(^)-^ AW )-.. 

where 22(0, z; A) is the gauge link, u(p', s'),u(p, s) are the Dirac spinors and 
we use the notation ^y a z a = z. 



xp+yr I \ xp-(l-y) r 
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a) b) c) 

Fig. 1. a) Parton picture in terms of y-DDs; 6, c) FM-type contributions. 

The parameter x in this representation has an evident interpretation: it 
characterizes which the fraction of the initial hadron momentum is carried 
by the active parton. 

In the nonforward case, we can use the double Fourier representation 
with respect to both (pz) and (rz): 

(p', s' | M0)zE(0, z; A)^ a (z) \p,s)\ 2 2 =0 (2) 
= j\y e -«(p*)-*y(«0 ! [ u(p',s')zu(p,s)F a (x,y;t) 

+ u(p',s') 4M u(p,s)K a (x,y;t) J 9{0 < x + y < 1) dx 

(p',s')u(p,s) Pe-^Vaiytfdy, 
Jo 



(zr) 
m 



where M is the nucleon mass and s, s' specify the nucleon polarization. We 
use the "hat" (rather than "slash") convention z = z^j^. The parametriza- 
tion of nonforward matrix elements must include both the nonflip term de- 
scribed here by the functions F a (x, y; t) and the spin-flip term characterized 
by the functions K a (x, y; t). 
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The parameters x, y tell us that the active parton carries the fractions x 
of the initial momentum p and the fraction y of the momentum transfer r. 
Using the approach [19] based on the a-representation analysis it is possible 
to prove [11] that double distributions F(x,y) have a natural property that 
both x and y satisfy the "parton" constraints < x < 1, < y < 1 for 
any Feynman diagram contributing to F(x,y). A less obvious restriction 
< x + y < 1 guarantees that the argument X = x + y( of the skewed 
distribution J-^(X) also changes between the limits The support area for 
the double distribution F a (x,y;t) is shown on Fig.2a. a 





a) 

Fig. 2. a) Support region and symmetry line y 
support region for a-DDs f(x, a). 



b) 

x/2 for y-DDs F(x,y;t); b) 



In principle, we cannot exclude also the possibility that the functions 
F(x, y; t) have singular terms at x = proportional to 5(x) or its deriva- 
tive^). Such terms have no projection onto the usual parton densities. We 
will denote them by Fm(x, y; t) — they may be interpreted as coming from 
the t-channel meson-exchange type contributions (see Fig. lb). In this case, 
the partons just share the plus component of the momentum transfer r: 
information about the magnitude of the initial hadron momentum is lost if 
the exchanged particle can be described by a pole propagator ~ \/{t — m 2 M ). 
Hence, the meson-exchange contributions to a double distribution may look 
like 

F+(x, y] t)~5(x)^M^7 or F M (x, y; t) ~ 8'(x) , e tc. , (3) 

% t m M t 

where <Pm(u) are the functions related to the distribution amplitudes of the 
relevant mesons M^. The two examples above correspond to x-even and 
x-odd parts of the double distribution F(x,y;t). Another type of terms in 
which the dependence on (pz) is lost can be produced by diagrams contain- 
ing a quartic pion vertex (Fig.lc). As shown by Polyakov and Weiss [20], 
such terms correspond to an independent (rz)u(j)', s')u(p, s)$((rz)) type 
contribution which can be parametrized by a single integral over y involv- 
ing an effective distribution amplitude ^(y;t). The meson-exchange terms 
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in F(x, y; t) and K(x, y; t) as well as Polyakov- Weiss terms are invisible in 
the forward limit, hence the existing knowledge of the usual parton densities 
cannot be used to constrain these terms. Later, describing the models for 
skewed distributions, we discuss only the "forward visible parts" of SPDs 
which are obtained by scanning the parts of the relevant DDs. 

Comparing Eq. (1) with the r = limit of the DD definition (2) gives 
the "reduction formulas" relating the double distribution F a (x,y;t = 0) to 
the quark and antiquark parton densities 



F a (x,y;t = 0)\ x>0 dy = f a (x) ; / F a (x,y;t = 0)\ x<0 dy = -fa(-x) . 



Hence, the positive- x and negative-x components of the double distribu- 
tion F a (x,y;t) can be treated as nonforward generalizations of quark and 
antiquark densities, respectively. If we define the "untilded" DDs by 

F a (x,y;t) = F a (x,y;t)\ x>0 ; F & (x, y; t) = -F a (-x, 1 - y; t)\ x<0 , (5) 

then x is always positive and the reduction formulas have the same form 



in both cases. The new antiquark distributions also "live" on the triangle 
< x,y < 1, < x + y < 1. Taking z in the lightcone "minus" direction, 
we arrive at the parton interpretation of functions F a ^{x,y\t) as probabil- 
ity amplitudes for an outgoing parton to carry the fractions xp + and yr + 
of the external momenta r and p. The double distributions F(x, y; t) are 
universal functions describing the flux of p + and r + independently of the 
ratio r + /p + . Note, that extraction of two separate components F a (x, y;t) 
and Fa(x, y; t) from the quark DD F a (x, y; t) as its positive-x and negative-x 
parts is unambiguous. 

Taking the O(z) term of the Taylor expansion gives the sum rules 





(4) 




(6) 
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relating the double distributions F a (x,y;t), K a (x,y;t) to the a-flavor com- 
ponents of the Dirac and Pauli form factors: 

Y,e a *f(t) = F 1 (t) , ]>> a F 2 a (t) = F 2 (i), (9) 

a a 

respectively. 

A common element of the reduction formulas given above is an integra- 
tion over y. Hence, it is convenient to introduce intermediate functions 

T a {x-t)= f 1 X F a (x,y;t)dy ; JC a (x;t) = C * K a (x,y;t) dy . (10) 

JO JO 

They satisfy the reduction formulas 

F a (x;t = 0) = f a (x) ; £ e a [' [F a (x; t) - T a {x- 1)\ dx = F 1 (t) (11) 

J2e a C [K a {x-t)-K^{x;t)\ dx = F 2 (t), (12) 

which show that these functions are the simplest hybrids of the usual parton 
densities and form factors. For this reason, one can call them nonforward 
parton densities (NDs) [21]. 

The spin-flip terms disappear only if r = 0. In the weaker r 2 = t = 
limit, they survive, e.g., F^O) = n a is the a-flavor contribution to the 
nucleon anomalous magnetic moment. In other words, the t = limit of 
the "magnetic" NDs exists: K, a (x;t = 0) = k a (x), and the integral 

Y]e a / [k a (x) - ka{x)} dx = k p (13) 

gives the anomalous magnetic moment of the proton. The knowledge of the 
x-moment of /c a (x)'s is needed to determine the contribution of the quark 
orbital angular momentum to the proton spin [8]. Since the i^T-type DDs are 
always accompanied by the = — p'^ factor, they are invisible in deep 
inelastic scattering and other inclusive processes related to strictly forward 
r = matrix elements. 

There are also parton-helicity sensitive double distributions G a (x,y;t) 
and P a (x,y;t). The first one reduces to the usual spin-dependent densities 
Af a (x) in the r = limit and gives the axial form factor F^t) after the 
x, y-integration. The second one involves hadron helicity flip and is related 
to the pseudoscalar form factor Fp(t). 
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xP+(l+a) r/2 / \ xP-(l- a) r/2 



P+r/2 P-r/2 



d) 

Fig. 3. Parton picture in terms of a-DDs. 

It is worth mentioning here that for a massive target (nuclcons in our 
case) there is a kinematic restriction [16] 

-t > Cm 2 H. (14) 

Hence, for fixed the formal limit i — > is not physically reachable. How- 
ever, many results (evolution equations being the most important example) 
obtained in the formal t = 0, M = limit are still applicable. 

To make the description more symmetric with respect to the initial and 
final hadron momenta, we can treat nonforward matrix elements as func- 
tions of (Pz) and (rz), where P = (p + p')/2 is the average hadron momen- 
tum. The relevant double distributions f a (x, a ; t) [which we will call a-DDs 
to distinguish them from y-DDs F(x, y; t)] are defined by 

(p'\M-z/2)zMz/2)\p) (15) 

= u(p')zu(p) f 1 dx f 1 e- ix(Pz) - iairz)/2 f a (x,a;t) da + 0(r) terms . 
J-i i-i+l^l 

The support area for f a (x,a;t) is shown in Fig.2b. Again, the usual for- 
ward densities f a (x) and fa{x) are given by integrating f a (x,a ; t = 0) over 
vertical lines x = const for x > and x < 0, respectively. Hence, we can 
split f a (x,a; t) into three components 

f a (x, a;t) = f a (x, a ; t) 9(x > 0) - fa{-x, -a ; t) 9(x < 0) + f M (x, a;t), 

(16) 

where fM(x,a; t) is a singular term with support at x = only. Due to 
hermiticity and time-reversal invariance properties of nonforward matrix 
elements, the a-DDs are even functions of a: 



f a (x,a;t) = f a (x,-a;t) 
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For our original y-DDs F a ^a(x,y;t), this corresponds to symmetry with re- 
spect to the interchange y <-> 1 — x — y ("Munich" symmetry, established 
in Ref. [22]). In particular, the functions (p^(y) for singular contributions 
F^(x,y;t) are symmetric <p^{y) = ^(1 — y) both for x-even and x-odd 
parts. The a-quark contribution into the flavor-singlet operator 

O s a (-z/2,z/2) = % -[^-z/2)zE{-z/2,z/2;A)^ a {z/2) - {z ^ -z}\ 

can be parametrized either by y-DDs F^(x,y;t) or by a-DDs (x, a; t) 
(p',s f \0'!(-z/2,z/2) \p,s)\ z2=0 (17) 

= u(p', S ')zu(p,s) ( dx [ I L-ix(pz)-i(y-l/2){rz) 

Jo Jo 2 V 

_ e ix(pz)+i(y-i/2)(rz)^ pS ( x ^ y;t ) dy + (r) terms 

= u(p',s')zu{p,s) f 1 dx f 1 e~ ix(Pz) ~ ia{rz)/2 f^{x,a;t)da + 0(r) . 
J-l J-l+\x\ 

In the second and third lines here we have used the fact that positive- x and 
negative-x parts in this case are described by the same untilded function 

Fa(x,y,t)\ XJ L = F a (x,y;t) + F a {x,y;t). 

The a-DDs /<f (x, a ; t) are even functions of a and, according to Eq. (17), 
odd functions of x: 

fa(x,a;t) = {f a (\x\,\a\;t) + f a (\x\,\a\;t)}s[gn(x) + f'^(x,a;t) . (18) 
Finally, the valence quark functions f^(x,a ; t) related to the operators 

OU-z/2,z/2) = ^[M-z/2)zE(-z/2,z/2;A^ a (z/2) + {z^-z}} 

are even functions of both a and x: 

f^(x,a;t) = f a (\x\, \a\;t) - fa(\x\, \a\;t) + f^(x,a;t) . (19) 
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3. Models for double and skewed distributions 

The reduction formulas and interpretation of the x-variable as the frac- 
tion of p (or P) momentum suggest that the profile of F(x,y) (or f(x,a)) 
in x-direction is basically determined by the shape of f(x). On the other 
hand, the profile in y (or a) direction characterizes the spread of the parton 
momentum induced by the momentum transfer r. In particular, since the 
a-DDs f(x,a) are even functions of a, it make sense to write 

f(x,a)=h(x,a)f(x), (20) 

where h(x, a) is an even function of a normalized by 

rl-|x| 

h(x,a)da = 1. (21) 

l+\x\ 



We may expect that the a-profile of h(x, a) is similar to that of a symmet- 
ric distribution amplitude (DA) ip(a). Since \a\ < 1 — \x\, to get a more 
complete analogy with DA's, it makes sense to rescale a as a = (1 — \x\)(3 
introducing the variable (3 with x-independent limits: — 1 < (3 < 1. The 
simplest model is to assume that the profile in the /3-direction is a universal 
function g(/3) for all x. Possible simple choices for g(/3) may be 5(/3) (no 
spread in /3-direction) , |(1 — /3 2 ) (characteristic shape for asymptotic limit 
of nonsinglet quark distribution amplitudes), j|(l — /? 2 ) 2 (asymptotic shape 
of gluon distribution amplitudes), etc. In the variables x,a, this gives 

*<«»(*, a) = 5(a) , hV{x,a) = ^ > 

r(2)/ . 15 [(1- |x|) 2 -a 2 ] 2 , . 

h(2) ^ = V, (i- N)5 • ^ 

These models can be treated as specific cases of the general profile function 
a) - H26 + 2) [(l-N) 2 - a 2 f 

l^«J - 2 2fe+i r 2 (6+1 ) (l_| x |)26+l ' W 

whose width is governed by the parameter b. 

The coefficient of proportionality ( = r + /p + (or £ = r + /2P + ) between 
the plus components of the momentum transfer and initial (or average) mo- 
mentum specifies the skewedness of the matrix elements. The characteristic 
feature implied by representations for double distributions [see, e.g., Eq.(2)] 
is the absence of the C-dependence in the DDs F(x, y) and ^-dependence in 
f(x,a). An alternative way to parametrize nonforward matrix elements of 
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a 



y 
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-1 




-1 



X 



a) 



Fig. 4. Integration lines for integrals relating SPDs and DDs. 



light-cone operators is to use £ (or £) and the total momentum fractions 
X = x + y( (or x = x + £a) as independent variables. Integrating each par- 
ticular double distribution over y gives the nonforward parton distributions 



where £ = 1 — C- The two components of NFPDs correspond to positive 
(X > and negative (X < Q values of the fraction X' = X — Q associated 
with the "returning" parton. As explained in refs. [9, 11], the second 
component can be interpreted as the probability amplitude for the initial 
hadron with momentum p to split into the final hadron with momentum 
(1 — ()p and the two-parton state with total momentum r = (p shared by 
the partons in fractions Yr and (1 — Y)r, where Y = X/Q. 

The relation between "untilded" NFPDs and DDs can be illustrated on 
the "DD-life triangle" defined by < x,y,x+y < 1 (seeFig.4a). Specifically, 
to get JT^(A), one should integrate F(x,y) over y along a straight line 
x = X — (y. Fixing some value of one deals with a set of parallel lines 
intersecting the x-axis at x = X. The upper limit of the y-integration is 
determined by intersection of this line either with the line x + y = 1 (this 
happens if X > Q) or with the y-axis (if X < (). The line corresponding to 
X = ( separates the triangle into two parts generating the two components 
of the nonforward parton distribution. 




0(X>() Fi(X-yC,y)dy + e(X<0 F t (X - y(, y) dy , 

Jo Jo 
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In a similar way, we can write the relation between OFPDs H(x, £; t) 
and the a-DDs f(x, a; t) 

/■i /-i-W 

H(x,£;t)= / cix / 5(x + £a — x) /(x, a; i) da . (25) 

J-l J-l+\x\ 

The delta-function in Eq.(25) specifies the line of integration in the {x,a} 
plane. For definiteness, we will assume below that £ is positive. 

Information contained in SPDs originates from two physically different 
sources: meson-exchange type contributions (X) coming from the sin- 
gular x = parts of DDs and the functions iF£(X), T^(X) obtained by 
scanning the i / parts of DDs F a (x,y), F a (x,y). The support of ex- 
change contributions is restricted to < X < £. Up to rescaling, the 
function T^ 1 (X) has the same shape for all £. For any nonvanishing X, 
these exchange terms become invisible in the forward limit ( — > 0. On the 
other hand, the support of functions J-?(X), THX) in general covers the 
whole < X < 1 region. Furthermore, the forward limit of such SPDs as 
J^' a (X) is rather well known from inclusive measurements. Hence, infor- 
mation contained in the usual (forward) densities f a (x), f a (x) can be used 
to restrict the models for JF?(X), T^(X). 

Let us consider SPDs constructed using simple models of DDs speci- 
fied above. In particular, the model f^°°\x,a) = 5(a) f(x) (equivalent to 
F(°°\x,y) = 5(y — x/2)/(x)), gives the simplest model H(°°>(x,£;t = 0) = 
f(x) in which OFPDs at t = have no ^-dependence. For NFPDs this gives 

i.e., NFPDs for non-zero ( are obtained from the forward distribution f(X) = 
J~<;=o(X) by shift and rescaling. 

In case of the 6=1 and 6 = 2 models, simple analytic results can be 
obtained only for some explicit forms of f(x). For the "valence quark" - 
oriented ansatz f^\x,a), the following choice of a normalized distribution 

/(i)( s ) = r(5 ~ a \ x -(l - x f (27) 
v ' 6T(1 - a) y ' v ' 

is both close to phenomenological quark distributions and produces a simple 
expression for the double distribution since the denominator (1 — x) 3 factor 
in Eq. (22) is canceled. As a result, the integral in Eq. (25) is easily 
performed and we get 



*v(*>0l|*l>€ = ^(l-f){ 



{2 - a )ai-x)(xl~ a + x 2 2 - a ) 



+ (f - £)( x 2-« _ x 2- 0) m + ( ~ ^ (28) 
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for \x\ > £ and 

H v&> OI|x|<€ = ^3 ( l ~ l) {*i~ a [( 2 " " *) + ^ " *)] 

+(£->-£) | (29) 

in the middle — £ < x < £ region. We use here the notation x± = (x + 
£)/(l + and X2 = (x — £)/(l — £) [23]. To extend these expressions onto 
negative values of £, one should substitute £ by |£|. One can check, however, 
that no odd powers of |£| would appear in the x N moments of H iv (x,£). 
Furthermore, these expressions are explicitly non-analytic for x = ±£. This 
is true even if a is integer. Discontinuity at x = ±£, however, appears only 
in the second derivative of H 1V (x,£), i.e., the model curves for H iv (x,£) 
look very smooth (see Fig. 5). The explicit expressions for NFPDs in this 
model were given in ref.[17]. The relevant curves are also shown in Fig.5. 




0.2 0.4 0.6 0.8 1 -1 -0.5 0.5 1 



Fig.5. Valence quark distributions: untilded NFPDs F«(x) (left) and OFPDs 
Hy(x,^) (right) with a = 0.5 for several values of (, namely, 0.1, 0.2, 0.4, 0.6, 
0.8 and corresponding values of £ = £/(2 — (). Lower curves correspond to larger 
values of (. 

For a = 0, the x > £ part of OFPD has the same x-dependence as its 
forward limit, differing from it by an overall ^-dependent factor only: 

H lv (x,O\a=o = ^ 1 ~_ l ^ 0(\x\>O + 2 * + ^~y /€ fl(|a|<fl- (30) 

The (1 — |x|) 3 behaviour can be trivially continued into the |x| < £ region. 
However, the actual behaviour of H iv (x, £)| a =o i n this region is given by 
a different function. In other words, H iv (x, £)| a =o can be represented as 
a sum of a function analytic at border points and a contribution whose 
support is restricted by |x| < £. It should be emphasized that despite its 
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DA-like appearance, this contribution should not be treated as an exchange- 
type term. It is generated by regular i/O part of DD, and, unlike (p(x/£)/£ 
functions changes its shape with £ and becomes very small for small £. 




Fig. 6. Singlet quark distribution Hg(x,£) for several £ values 0.1, 0.25, 0.4. 
Lower curves correspond to larger values of £. Forward distribution is modeled 
by (1 - x) 3 /x. 

For the singlet quark distribution, the a-DDs f s (x,a) should be odd 
functions of x. Still, we can use the model like (27) for the x > part, 
but take f S (x,a)\ x ^Q = Af^\\x\, a) sign(x). Note, that the integral (25) 
producing H (x,£) in the \x\ < £ region would diverge for a — > x/£ if 
a > 1, which is the usual case for standard parametrizations of singlet quark 
distributions for sufficiently large Q 2 . However, due to the antisymmetry 



of f (x, a) wrt x — ► — x and its symmetry wrt a 



-a 



the singularity at 



a = x/£ can be integrated using the principal value prescription which in 
this case produces the x — ► — x antisymmetric version of Eqs.(28) and (29). 
For a = 0, its middle part reduces to 



H lS (x,£)\\x\<t,a=0 = 2x 



3£ 2 - 2x 2 i - x 2 



The shape of singlet SPDs in this model is shown in Fig. 6 



(31) 



4. SPDs and deeply virtual Compton scattering 

In the lowest order, the DVCS amplitude T^ u (p,q,q') is given by two 
handbag diagrams. In particular, the invariant amplitude containing the T 
functions is given by 



T F (p,q,q') =^ e a / 

n JO 



_X -C + ie + X-ie 



T* + ~ a (X;t)dX . (32) 
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An important feature of the DVCS amplitude is that for large Q 2 and fixed 
t it depends only on the ratio Q 2 /2(pq) = XBj = (■ DVCS is an exclusive 
process exhibiting the Bjorken scaling. Note that the imaginary part of 
the DVCS amplitude is proportional to J-^ +a ((,', t) . In this function, the 
parameter £ appears twice: first as the skewedness of the process and then 
as the fraction X = ( at which the imaginary part is generated. 

One may ask which Q 2 are large enough to ensure the dominance of the 
lowest-twist handbag contribution. In DIS, approximate Bjorken scaling 
starts at Q 2 ~ 2 GeV 2 . Another example is given by the exclusive process 
7(<?i)7*te) - ► 71-0 studied on e + e~ colliders. If one of the photons is highly 
virtual q 2 = —Q 2 while another is (almost) real q\ ~ 0, the process is kine- 
matically similar to DVCS. In the leading order, the F^* n o(Q 2 ) transition 
form factor is given by a handbag diagram again. The recent measurements 
by CLEO [24] show that the pQCD prediction F 77 * 7r o(Q 2 ) ~ 1/Q 2 again 
works starting from Q 2 ~ 2 GeV 2 . The 77*71"° vertex (for a virtual pion) 
can be also measured on a fixed-target machine like CEBAF in which case 
it is just a part of the DVCS amplitude corresponding to the 4th skewed 
distribution V^(X, t) (which is related to the pseudoscalar form factor Gp(t) 
of the nucleon). Hence, CLEO data give an evidence that DVCS may be 
handbag-dominated for Q 2 as low as 2 GeV 2 . 

The main problem for studying DVCS is the contamination by the 
Bethe-Heitler process in which the final photon is emitted from the ini- 
tial or final electron. The Bethe-Heitler amplitude is enhanced at small t. 
On the other hand, the virtual photon flux for fixed Q 2 and XBj increases 
when the electron beam energy increases. Hence, the energy upgrade would 
make the DVCS studies at Jefferson Lab more feasible. Experimental as- 
pects of virtual Compton scattering studies at Jefferson Lab are discussed 
in Refs. [25]. 

The skewed parton distributions can be also measured in hard meson 
electroproduction processes [11, 12, 22, 25]. The leading-twist pQCD contri- 
bution in this case involves a one-gluon exchange, which means that the hard 
subprocess is suppressed by a s /ir ~ 0.1 factor. The competing soft mecha- 
nism corresponds to a triple overlap of hadronic wave functions and has a 
relative suppression M 2 /Q 2 by a power of Q 2 , with M 2 ~ 1 GeV 2 being a 
characteristic hadronic scale. Hence, to clearly see the one-gluon-exchange 
signal one needs Q 2 above 10 GeV 2 . Numerical pQCD-based estimates and 
comparison of DVCS and hard meson electroprodution cross sections can 
be found in Ref. [25]. 
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5. SPD enhancement factor 

The imaginary part of hard exclusive meson electroproduction amplitude 
is determined by the skewed distributions at the border point. For this 
reason, the magnitude of ^(C) [or -ff (£,£)] and its relation to the forward 
densities f(x) has a practical interest. This example also gives a possibility 
to study the sensititivity of the results to the choice of the profile function. 
Assuming the infinitely narrow weight p(a) = 6(a), we have ^(X) = f(X— 
C/2) + ... and H{x,^) = f(x). Hence, both ^(C) and H(^,^) are given 
by f(xsj/2) since C = xb-j and £ = XBj/2 + .... Since the argument of 
f(x) is twice smaller than in deep inelastic scattering, this results in an 
enhancement factor. In particular, if f(x) ~ x~ a for small x, the ratio 
TZ(C) = ^(C)//(C) is 2 a . The use of a wider profile function p{a) produces 
further enhancement. For example, taking the normalized profile function 

and /(x) ~ x~ a we get 

-.(&) 



v(b)(n = jy^l = r(2fc + 2)r(b-a + i) 

/(C) r(26-a + 2)r(6 + l) 



(34) 



which is larger than 2 a for any finite b and < a < 2. The 2 a enhancement 
appears as the b — > oo limit of Eq.(33). For small integer b, Eq.(33) reduces 
to simple formulas obtained in refs. [17, 18]. For b = 1, we have 



■^f =1) (C) l 



/(C) (l-a/2)(l-a/3) 



(35) 



which gives the factor of 3 for the enhancement if a = 1. For 6 = 2, the 
ratio (33) becomes 



^ b=2 \0 i 



/(C) (l-a/3)(l-a/4)(l-a/5)' 



(36) 



producing a smaller enhancement factor 5/2 for a = 1. Calculating the 
enhancement factors, one should remember that the gluon SPD T^{X) re- 
duces to Xfg(X) in the C = limit. Hence, to get the enhancement factor 
corresponding to the f g (x) ~ x~ A small-x behavior of the forward gluon 
density, one should take a = A — 1 in Eq.(33), i.e., despite the fact that 
the 1/x behavior of the singlet quark distribution gives the factor of 3 for 
the ^(^(C) ratio, the same shape of the gluon distribution results in no 
enhancement. 
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Due to evolution, the effective parameter a characterizing the small-x 
behavior of the forward distribution is an increasing function of Q 2 . As a 
result, for fixed b, the TZ^(Q ratio increases with Q 2 . In general, the profile 
of f(x,a) in the a-direction is also affected by the pQCD evolution. In 
particular, in ref. [17] it was shown that if one takes an ansatz corresponding 
to an extremely asymmetric profile function p(a) ~ <5(l + a), the shift of the 
profile function to a more symmetric shape is clearly visible in the evolution 
of the relevant SPD. Recently, it was demonstrated [26, 27] that evolution to 
sufficiently large Q 2 enforces a direct relation b = a between the parameter 
a characterizing the small-x behavior of DDs and the parameter b governing 
the shape of their a profile. This gives 

n (b=a) (n r ( 2a + 2 ) r 3 7 <| 

K Kj "r(a + 2)r(a + i) [67) 

for the TZ{() ratio. For a = 1, e.g., the SPD enhancement factor in this case 
equals 3. 



6. Compton scattering 

6.1. General Compton amplitude 

The Compton scattering in its various versions provides a unique tool for 
studying hadronic structure. The Compton amplitude probes the hadrons 
through a coupling of two electromagnetic currents and in this aspect it can 
be considered as a generalization of hadronic form factors. In QCD, the 
photons interact with the quarks of a hadron through a vertex which, in 
the lowest approximation, has a pointlike structure. However, in the soft 
regime, strong interactions produce large corrections uncalculable within 
the perturbative QCD framework. To take advantage of the basic pointlike 
structure of the photon-quark coupling and the asymptotic freedom feature 
of QCD, one should choose a specific kinematics in which the behavior of 
the relevant amplitude is dominated by short (or, being more precise, light- 
like) distances. The general feature of all such types of kinematics is the 
presence of a large momentum transfer. For Compton amplitudes, there 
are several situations when large momentum transfer induces dominance of 
configurations involving lightlike distances: 

i) both photons are far off-shell and have equal spacelike virtuality: virtual 
forward Compton amplitude, its imaginary part determines structure func- 
tions of deep inelastic scattering (DIS); 

ii) initial photon is highly virtual, the final one is real and the momentum 
transfer to the hadron is small: deeply virtual Compton scattering (DVCS) 
amplitude; 
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Hi) both photons are real but the momentum transfer is large: wide-angle 
Compton scattering (WACS) amplitude. 

Our main statement made in ref. [21] is that, at accessible momentum 
transfers \t\ < 10 GeV 2 , the WACS amplitude is dominated by handbag dia- 
grams, just like in DIS and DVCS. In the most general case, the nonpertur- 
bative part of the handbag contribution is described by double distributions 
(DDs) F(x,y;t),G(x,y;t), etc., which can be related to the usual parton 
densities f(x), A/(x) and nucleon form factors like Fi (i) , (i) ■ Among 
the arguments of DDs, x is the fraction of the initial hadron momentum 
carried by the active parton and y is the fraction of the momentum transfer 
r. The description of WACS amplitude simplifies when one can neglect the 
y-dependence of the hard part and integrate out the y-dependence of the 
double distributions. In that case, the long-distance dynamics is described 
by nonforward parton densities (NDs) J-(x;t),Q(x;t), etc. The latter can 
be interpreted as the usual parton densities f(x) supplemented by a form 
factor type t-dependence. We proposed in [21] a simple model for the rele- 
vant NDs which both satisfies the relation between J-{x\ t) and usual parton 
densities f(x) and produces a good description of the i*i(t) form factor up to 
t ~ —10 GeV 2 . We have used this model to calculate the WACS amplitude 
and obtained the results which are rather close to existing data. 

6.2. Deep inelastic scattering 

The forward virtual Compton amplitude whose imaginary part gives 
structure functions of deep inelastic scattering (see, e.g., [1]) is the classic 
example of a light cone dominated Compton amplitude. In this case, the 
"final" photon has momentum q' = q coinciding with that of the initial 
one. The momenta p, p' of the initial and final hadrons also coincide. The 
total cm energy of the photon-hadron system s = (p + q) 2 should be above 
resonance region, and the Bjorken ratio xbj = Q 2 /2(pq) is finite. The light 
cone dominance is secured by high virtuality of the photons: — q 2 = Q 2 > 1 
GeV 2 . In the large-Q 2 limit, the leading contribution in the lowest a s order 
is given by handbag diagrams in which the perturbatively calculable hard 
quark propagator is convoluted with parton distribution functions f a (x) 
(a = u,d,s, . . .) which describe/parametrize nonperturbative information 
about hadronic structure. 

6.3. Deeply virtual Compton scattering 

The condition that both photons are highly virtual may be relaxed by 
taking a real photon in the final state. Keeping the momentum transfer 
t = (p — p') 2 to the hadron as small as possible, one arrives at kinematics of 
the deeply virtual Compton scattering (DVCS) the importance of which was 
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recently emphasized by X. Ji [8] (see also [9] ) . Having large virtuality Q 2 of 
the initial photon is sufficient to guarantee that in the Bjorken limit the lead- 
ing power contributions in 1/Q 2 are generated by the strongest light cone 
singularities [8, 11, 28, 29], with the handbag diagrams being the starting 
point of the a s expansion. The most important contribution to the DVCS 
amplitude is given by a convolution of a hard quark propagator and a non- 
perturbative function describing long-distance dynamics, which in the most 
general case is given by double distributions (DDs) F(x, y; t), G(x, y;t), . . . 
I']. 

In the DVCS kinematics, \t\ is assumed to be small compared to Q 2 , and 
for this reason the t- and m^-dependence of the short-distance amplitude in 
refs. [8, 9, 11, 23] was neglected 1 . This is equivalent to approximating the 
active parton momentum k by its plus component alone: k — > xp + + yr + . 

7. Modeling NDs 

Our final goal in the present paper is to get an estimate of the hand- 
bag contributions for the large-t real Compton scattering. Since the initial 
photon in that case is also real: Q 2 = (and hence xb-j = 0), it is nat- 
ural to expect that the nonperturbative functions which appear in WACS 
correspond to the ( = limit of the skewed parton distributions 2 ^(x^t). 
It is easy to see from Eq.(10) that in this limit the SPDs reduce to the 
nonforward parton densities F a (x;t) introduced above: 

^= (x;t) = F a {x-t) . (38) 

1 One should not think that such a dependence is necessarily a higher twist effect: the 
lowest twist contribution has a calculable dependence on t and nip analogous to the 
Nachtmann-Georgi-Politzer 0(irip/Q 2 ) target mass corrections in DIS [30, 31]. 

2 Provided that one can neglect the t-dependence of the hard part. 



Skewed Parton Distributions 



3665 



Note that NDs depend on "only two" variables x and t, with this dependence 
constrained by reduction formulas (11), (12). Furthermore, it is possible to 
give an interpretation of nonforward densities in terms of the light-cone 
wave functions. 

Consider for simplicity a two-body bound state whose lowest Fock com- 
ponent is described by a light cone wave function ^(x, k±). Choosing a 
frame where the momentum transfer r is purely transverse r = rj_, we can 
write the two-body contribution into the form factor as [32] 

F^(t) = J\x J *%x,k ± +xr ± )*(x,k ± )^±, (39) 
where x = 1 — x. Comparing this expression with the reduction formula 




a) b) 



Fig. 8. a) Structure of the effective two-body contribution to form factor in the light 
cone formalism, b) Form factor as an x-integral of nonforward parton densities. 

(11), we conclude that 

fW(x,t)= J y*( x ,k ± + xr ± )V(x,k ± )^± (40) 

is the two-body contribution into the nonforward parton density .F(x,i). 
Assuming a Gaussian dependence on the transverse momentum k± (cf. [32]) 

V(x,k ± ) = §{x)e- k l/ 2xix2 , (41) 

we get 

J* tb )(x,t) = f( tb \x)e st ^ xX2 , (42) 

where 

f( tb \x) = ^&(x)=FW(x,t = 0) (43) 
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is the two-body part of the relevant parton density. Within the light-cone 
approach, to get the total result for either usual f(x) or nonforward par- 
ton densities T{x,t), one should add the contributions due to higher Fock 
components. By no means these contributions are small, e.g., the valence 
du contribution into the normalization of the tt + form factor at t = is less 
than 25% [32]. In the absence of a formalism providing explicit expressions 
for an infinite tower of light-cone wave functions we choose to treat Eq.(42) 
as a guide for fixing interplay between the t and x dependences of NDs and 
propose to model them by 

F a ( X ,t) = f a (x)e St ^ xX2 = I' e -{klHk±+xr^)/2x^ d 2 k± _ {u) 

irxx\ z J 

The functions f a (x) here are the usual parton densities assumed to be taken 
from existing parametrizations like GRV, MRS, CTEQ, etc. In the t = 
limit (recall that t is negative) this model, by construction, satisfies the 
first of reduction formulas (11). Within the Gaussian ansatz (44), the basic 
scale A specifies the average transverse momentum carried by the quarks. 
In particular, for valence quarks 

2 \a _ ^ I " „- rval 



J^a Jo 

where N u = 2, Nd = 1 are the numbers of the valence a-quarks in the proton. 

To fix the magnitude of A, we use the second reduction formula in (11) 
relating J- a (x, t)'s to the F\ (t) form factor. To this end, we take the following 
simple expressions for the valence distributions 

fl a \x) = 1.89x _a4 (l - x) 3 - 5 (l + 6x) , (46) 

f™\ x ) = 0.54x- a6 (l - x) 4 - 2 (l + 8s) . (47) 

They closely reproduce the relevant curves given by the GRV parametriza- 
tion [33] at a low normalization point Q 2 ~ 1 GeV 2 . The best agreement 
between our model 

F? oit (t) = C k fl a \x) + e d fr\x)} e^ AxX2 dx (48) 
Jo 1 J 

and experimental data [34] in the moderately large t region 1 GeV 2 < \t\ < 
10 GeV 2 is reached for A 2 = 0.7 GeV 2 (see Fig. 9). This value gives a 
reasonable magnitude 

(k 2 ± ) u = (290 MeV) 2 , {k\) d = (250 MeV) 2 (49) 
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Fig. 9. Ratio F[(t)/D(t) of the F[(t) form factor to the dipole fit D(t) = 1/(1 - 
i/0.71GeV 2 ) 2 . Curve is based on Eq. (47) with A 2 = 0.7 GcV 2 . Experimental 
data are taken from ref. [34] . 



for the average transverse momentum of the valence u and d quarks in the 
proton. 

Similarly, building a model for the parton helicity sensitive NDs G a (x, t) 
one can take their t = shape from existing parametrizations for spin- 
dependent parton distributions A/ a (x) and then fix the relevant A param- 
eter by fitting the GU(i) form factor. The case of hadron spin-flip distri- 
butions JC a (x,t) and V a (x;t) is more complicated since the distributions 
k a (x), p a {x) are unknown. 

At t = 0, our model by construction gives a correct normalization (t = 
0) = 1 for the form factor. However, if one would try to find the derivative 
(d/dt)F^(t) at t = by expanding the exponential exp[xt/x\ 2 ] into the 
Taylor series under the integral (48), one would get a divergent expression. 
An analogous problem is well known in applications of QCD sum rules to 
form factors at small t [35, 36, 37, 38]. The divergence is related to the 
long-distance propagation of massless quarks in the i-channel. Formally, 
this is revealed by singularities starting at t = 0. However, Ff(t) should 
not have singularities for timelike t up to 4m: 2 ., with the p- meson peak 
at t = m 2 p ~ 0.6 GeV 2 being the most prominent feature of the t-channel 
spectrum. Technically, the singularities of the original expression are singled 
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out into the bilocal correlators [39] which are substituted by their realistic 
version with correct spectral properties (usually the simplest model with 
p and p' terms is used). An important point is that such a modification 
is needed only when one calculates form factors in the small-i region: for 
— t > 1 GeV 2 , the correction terms should vanish faster than any power of 
1/t [37]. In our case, the maximum deviation of the curve for Ff(t) given 
by Eq.(48) from the experimental data in the small-i region —t < 1 GeV 2 
is 15%. Hence, if one is willing to tolerate such an inaccuracy, one can use 
our model starting with t = 0. 

Our curve is within 5% from the data points [34] for 1 GeV 2 < — t < 6 
GeV 2 and does not deviate from them by more than 10% up to 9 GeV 2 . 
Modeling the i-dependence by a more complicated formula (e.g., assuming 
a slower decrease at large t, and/or choosing different A's for u and d quarks 
and/or splitting NDs into several components with different A's, etc., see 
Ref.[40] for an example of such an attempt) or changing the shape of parton 
densities f a (x) one can improve the quality of the fit and extend agreement 
with the data to higher t. Such a fine-tuning is not our goal here. We 
just want to emphasize that a reasonable description of the Fi(t) data in 
a wide region 1 GeV 2 < \t\ < 10 GeV 2 was obtained by fixing just a single 
parameter A reflecting the proton size. Moreover, we could fix A from the 
requirement that (k^_) ~ (300 MeV) 2 and present our curve for Fi(t) as a 
successful prediction of the model. We interpret this success as an evidence 
that the model correctly catches the gross features of the underlying physics. 

Since our model implies a Gaussian dependence on the transverse mo- 
mentum, it includes only what is usually referred to as an overlap of soft 
wave functions. It completely neglects effects due to hard pQCD gluon 
exchanges generating the power-law 0((a s /7r) 2 /i 2 ) tail of the nonforward 
densities at large t. It is worth pointing out here that though we take nonfor- 
ward densities J- a (x, t) with an exponential dependence on t, the form 
factor in our model has a power-law asymptotics Ff oit (t) ~ (— 4A 2 /i) n+1 dic- 
tated by the (1 — x) n behavior of the parton densities for x close to 1. This 
connection arises because the integral (48) over x is dominated at large t 
by the region x ~ 4A 2 /|i|. In other words, the large-t behavior of Fi(t) in 
our model is governed by the Feynman mechanism [1]. One should realize, 
however, that the relevant scale 4A 2 = 2.8 GeV 2 is rather large. For this 
reason, when \t\ < 10 GeV 2 , it is premature to rely on asymptotic estimates 
for the soft contribution. Indeed, with n = 3.5, the asymptotic estimate is 
~ t 4 5 , in an apparent contradiction with the ability of our curve 
to follow the dipole behavior. The resolution of this paradox is very sim- 
ple: the maxima of nonforward densities J- a (x,t) for |i| < 10 GeV 2 are at 
rather low x- values x < 0.5. Hence, the x-integrals producing Ff oft (t) are 
not dominated by the x ~ 1 region yet and the asymptotic estimates are 
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not applicable: the functional dependence of Ff (t) in our model is much 
more complicated than a simple power of 1/t. 

The fact that our model closely reproduces the experimentally observed 
dipole-like behavior of the proton form factor is a clear demonstration that 
such a behavior may have nothing to do with the quark counting rules 
Ff(t) ~ 1/t 2 [41, 42] valid for the asymptotic behavior of the hard gluon 
exchange contributions. Our explanation of the observed magnitude and 
the t-dependence of F± (t) by a purely soft contribution is in strong contrast 
with that of the hard pQCD approach to this problem. 

8. Wide-angle Compton scattering 

With both photons real, it is not sufficient to have large photon energy to 
ensure short-distance dominance: large-s, small-t region is strongly affected 
by Regge contributions. Hence, having large \t\ > 1 GeV 2 is a necessary 
condition for revealing short-distance dynamics. 

The simplest contributions for the WACS amplitude are given by the 
s- and u-channel handbag diagrams Fig. 7b, c. The nonperturbative part in 
this case is given by the proton DDs which determine the t-dependence of 
the total contribution. Just like in the form factor case, the contribution 
dominating in the formal asymptotic limit s, \t\, \u\ — » oo, is given by dia- 
grams corresponding to the pure SD regime, see Fig. 10a. The hard subraph 
then involves two hard gluon exchanges which results in a suppression fac- 
tor (a s /7r) 2 ~ 1/100 absent in the handbag term. The total contribution of 
all two-gluon exchnange contributions was calculated by Farrar and Zhang 
[43] and then recalculated by Kronfeld and Nizic [44]. A sufficiently large 
contribution is only obtained if one uses humpy DAs and 1/k 2 propagators 
with no finite-size effects included. Even with such propagators, the WACS 
amplitude calculated with the asymptotic DA is negligibly small [45] com- 
pared to existing data. As argued in ref.[21], the enhancements generated 
by the humpy DAs should not be taken at face value both for form fac- 
tors and wide-angle Compton scattering amplitudes. For these reasons, we 
ignore hard contributions to the WACS amplitude as negligibly small. 

Another type of configurations containing hard gluon exchange is shown 
in Fig. 10b. There are also the diagrams with photons coupled to differ- 
ent quarks ("cat's ears", Fig.lOc). Such contributions have both higher 
order and higher twist. This brings in the a s /ir factor and an extra 1/s 
suppression. The latter is partially compensated by a slower fall-off of the 
four-quark DDs with t since only one valence quark should change its mo- 
mentum. 

For simplicity, we neglect all the suppressed terms and deal only with the 
handbag contributions Fig.7b,c in which the highly virtual quark propagator 
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connecting the photon vertices is convoluted with proton DDs parametrizing 
the overlap of soft wave functions. Since the basic scale 4A 2 characterizing 
the t-dependence of DDs in our model is 2.8 GeV 2 , while existing data are all 
at momentum transfers t below 5 GeV 2 , we deal with the region where the 
asymptotic estimate (Feynman mechanism) for the overlap contribution is 
not working yet. In the coordinate representation, the sum of two handbag 
contributions to the Compton amplitude can be written as 



+M-z/2)YS c (-z)YM^/mP)d 4 z 



(50) 



where Q = (q + q')/2 and S°(z) = iz/2n 2 (z 2 ) 2 is the hard quark propagator 
(throughout, we use the "hat" notation z = z a ~f a ). The summation over 
the twist-0 longitudinal gluons adds the usual gauge link between the ip,ip 
fields which we do not write down explicitly (gauge link disappears, e.g., in 
the Fock-Schwinger gauge z a A a (z) = 0). Because of the symmetry of the 
problem, it is convenient to use P = (p + p')/2 (cf. [8]) and r = p — p' 
as the basic momenta. Applying the Fiertz transformation and introducing 
the double distributions by 



(p'\M-z/2)zMz/2)\p) =u{p')zu{p) / dx 



-i(kz) fa 



f a (x,a;t) 



_ e i(kz) fa 



f a (x,a;t) dy + u(p')(zf — fz)u(p) 



4m 



(51) 



dx 



e - i{kz) k a (x, a; t) - e i{kz) k a (x, a; t)] da + 0{z 2 ) terms 



(we use here the shorthand notation k = xP + ar/2) and similarly for the 
parton helicity sensitive operators 

(p'\M-z/2)z l5 Mz/2)\p) =u(p')h 5 u(p) f 1 dx f \e- l ^g a (x,a;t) 

JO J-x 1 
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+e l ^g~ a (x,a;t)} da + M u(p') lb u{p) (52) 
x r \e- i ^p a (x,y;t) + e i ^ P d (x,a;t)] da + 0{z 2 ) terms, 

JO J-S L J 

we arrive at a leading-twist QCD parton picture with a-DDs serving as func- 
tions describing long-distance dynamics. The a-DDs f a (x,a;t), etc., are 
related to the original y-DDs F a (x, y; t) by the shift y = (1 — x + a)/2. Inte- 
grating f(x,a;t) over a one obtains the same nonforward densities T{x\ t). 
The hard quark propagators for the s and u channel handbag diagrams in 
this picture look like 



and 



xP + ar/2 + Q xP + ar/2 + Q 

(xP + ar/2 + Q) 2 xs - (x 2 - a 2 )t/4 + x 2 m 2 

xP + ar/2-Q xP + af/2-Q 



(xP + ar/2 — Q) 2 xu — (x 2 — a 2 )t/A + x 2 m 2 



(53) 



(54) 



p 



respectively. We denote s = 2{pq) = s—m 2 and u = —2(pq r ) = u—m 2 . Since 
a-DDs are even functions of a [22], the af terms in the numerators can be 
dropped. It is legitimate to keep 0(m 2 ) and 0(t) terms in the denominators: 
the dependence of hard propagators on target parameters m 2 and t can be 
calculated exactly because of the effect analogous to the ^-scaling in DIS [31] 
(see also [46]). Note that the t-correction to hard propagators disappears 
in the large-t limit dominated by the x ~ 1 integration. The ^corrections 
are the largest for y = 0. At this value and for x = 1/2 and t = u (cm 
angle of 90°), the t-term in the denominator of the most important second 
propagator is only 1/8 of the u term. This ratio increases to 1/3 for x = 1/3. 
However, at nonzero a-values, the t-corrections are smaller. Hence, the t- 
corrections in the denominators of hard propagators can produce 10% — 20% 
effects and should be included in a complete analysis. Here, we consider an 
approximation in which these terms are neglected and hard propagators are 
given by y-independent expressions (xP + Q)/xs and {xP + Q)/xu. As a 
result, the a-integration acts only on the DDs f(x,a;t) and converts them 
into nonforward densities F(x,t). The latter would appear then through 
two types of integrals 

f 1 F a (x,t)dx = F?(t) and C F a {x,t) — = R°(t), (55) 

JO JO x 

and similarly for /C, Q, V . The functions Ff(t) are the flavor components of 
the usual Fi(t) form factor while Rf(t) are the flavor components of a new 
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form factor specific to the wide-angle Compton scattering. In the formal 
asymptotic limit \t\ — > oo, the x integrals for Ff(t) and Rf(t) are both 
dominated in our model by the x ~ 1 region: the large-i behavior of these 
functions is governed by the Feynman mechanism and their ratio tends to 
1 as \t\ increases (see Fig. 11a). However, due to large value of the effective 
scale 4 A 2 = 2.8 GeV 2 , the accessible momentum transfers t < 5 GeV 2 are 
very far from being asymptotic. 

In Fig. lib we plot J^&t) and F u (x;t)/x at t = -2.5 GeV 2 . It is 
clear that the relevant integrals are dominated by rather small x values 
x < 0.4 which results in a strong enhancement of Ri(t) compared to F"(t) 
for \t\ < 5 GeV 2 . Note also that the (p'\ . . . xP . . . \p) matrix elements 
can produce only t as a large variable while (p'| . . . Q . . . \p) gives s. As 
a result, the enhanced form factors Rf(t) are accompanied by extra s/t 
enhancement factors compared to the Ffft) terms. In the cross section, 
these enhancements are squared, i.e., the contributions due to the non- 
enhanced form factors (t) are always accompanied by t 2 /s 2 factors which 
are smaller than 1/4 for cm angles below 90°. Because of double suppression, 
we neglect Ff(t) terms in the present simplified approach. 




Fig. 11. a) Ratio R^(t)/F^(t); b) Functions F u (x;t) (solid line) and T u {x;t)jx 
(dashed line) at t = -2.5 GeV 2 . 

The contribution due to the JC functions appears through the flavor 
components F%{t) of the F-^it) form factor and their enhanced analogues 
R^if)- The major part of contributions due to the /C-type NDs appears 
in the combination Rj(t) - (t/ '4m 2 )i?|(t). Experimentally, F 2 (i)/Fi(t) « 
lGeV 2 /|t|. Since R2/F2 ~ R\/F 1 ~ l/(x), R 2 {t) is similarly suppressed 
compared to Ri(t), and we neglect contributions due to the R^if) form 
factors. We also neglect here the terms with another spin-flip distribution 
V related to the pseudoscalar form factor Gp(t) which is dominated by 
the t-channel pion exchange. Our calculations show that the contribution 
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due to the parton helicity sensitive densities Q a is suppressed by the factor 
t 2 /2s 2 compared to that due to the T a densities. This factor only reaches 
1/8 for the cm angle of 90°, and hence the Q a contributions are not very 
significant numerically. For simplicity, we approximate Q a (x,t) by J- a (x,t). 
After these approximations, the WACS cross section is given by the product 



da 2ira 2 



dt s 



Z.2 



(pq) + ivi') 



(jpq 1 ) (pq) 



R 2 (t) , (56) 



of the Klein-Nishina cross section (in which we dropped 0(m 2 ) and 0(m 4 ) 
terms) and the square of the Ri(t) form factor 

Ri(t) = J2^a[RUt) + Ri(t)] ■ (57) 

a 

In our model, R\(t) is given by 



Ri(t) = C 

Jo 



el f v u al (x)+e 2 f T \ x )+2{ei+ei+ei) f ea (x) 



(58) 



We included here the sea distributions assuming that they are all equal 
f sea (x) = fu £ ds( x ) = fuds( x ) an d using a simplified parametrization 

f sea (x) =0.5x~°- 75 (l-x) 7 (59) 

which accurately reproduces the GRV formula for Q 2 ~ 1 GeV 2 . Due to 
suppression of the small-x region by the exponential exp[xt/4xA 2 ], the sea 
quark contribution is rather small (~ 10%) even for — t ~ 1 GeV 2 and is 
invisible for -t > 3 GeV 2 . 

Comparison with existing data [47] is shown in Fig. 12. Our curves fol- 
low the data pattern but are systematically lower by a factor of 2, with 
disagreement becoming more pronounced as the scattering angle increases. 
Since we neglected several terms each capable of producing up to a 20% 
correction in the amplitude, we consider the agreement between our curves 
and the data as encouraging. The most important corrections which should 
be included in a more detailed investigation are the t-corrections in the de- 
nominators of hard propagators and contributions due to the "non-leading" 
/C, Q, V nonforward densities. The latter, as noted above, are usually ac- 
companied by t/s and t/u factors, i.e., their contribution becomes more 
significant at larger angles. The i-correction in the most important hard 
propagator term l/[xu — (x 2 — a 2 )t/A + x 2 m 2 ] also enhances the amplitude 
at large angles. 

The angular dependence of our results for the combination s 6 (da/dt) 
is shown on Fig.13. All the curves for initial photon ehergies 2,3,4,5 and 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 

-t (GeV 2 ) 



Fig. 12. WACS cross section versus t: comparison of results based on Eq.(55) with 
experimental data. 

6 GeV intersect each other at 9 cm ~ 60°. This is in good agreement with 
experimental data of ref.[47] where the differential cross section at fixed cm 
angles was fitted by powers of s: da/dt ~ s~ n ( e ) with n cxp (60°) = 5.9 ± 0.3. 
Our curves correspond to n soft (60°) « 6.1 and n soft (90°) w 6.7 which also 
agrees with the experimental result n exp (90°) = 7.1 ± 0.4. 

This can be compared with the scaling behavior of the asymptotic hard 
contribution: modulo logarithms contained in the a s factors, they have 
a universal angle-independent power n hard (#) = 6. For 6 cm = 105°, the 
experimental result based on just two data points is n exp (105°) = 6.2 ± 1.4, 
while our model gives n soft (105°) « 7.0. Clearly, better data are needed to 
draw any conclusions here. 

9. Conclusions 

The hard exclusive electroproduction processes provide new information 
about hadronic structure accumulated in skewed parton distributions. The 
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Fig. 13. Angular dependence of the combination s 6 (da/dt). 



SPDs are universal hybrid functions having the properties of parton densi- 
ties, hadronic form factors and distribution amplitudes. They give a unified 
description of various hard exclusive and inclusive reactions. The basic 
supplier of information about skewed parton distributions is deeply virtual 
Compton scattering which offers a remarkable example of Bjorken scaling 
phenomena in exclusive processes. Furthermore, wide-angle real Compton 
scattering is an ideal tool to test angle-dependent scaling laws characteristic 
for soft overlap mechanism. 
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Fig. 14. s-dependence of the combination s 6 da/dt for 9 = 60° (dotted line), 9 = 90° 
(dashed line) and 9 — 105° (solid line). 
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